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Using the noncommutative Kubo formula for aperiodic solids and a recently developed numerical 
implementation, we study the conductivity o and resistivity p tensors as functions of Fermi level Ep 
and temperature T for models of strongly disordered Chern insulators. The formalism enabled us to 
converge the transport coefficients at temperatures low enough to enter the quantum critical regime at 
the Chern-to-trivial insulator transition. We find that the p^x-curves at different temperatures intersect 
each other at one single critical point, and that they obey a single-parameter scaling law with an 
exponent close to the universally accepted value for the unitary symmetry class. However, when 
compared with the established experimental facts on the plateau-insulator transition in the Integer 
Quantum Hall Effect, we find a universal critical conductance cr^^ twice as large, an ellipse rather than 
a semi-circle law, and absence of the Quantized Hall Insulator phase. 

PACS numbers: 72.25.-b, 72.10.Fk, 73.20.Jc, 73.43.-f 



1. INTRODUCTION 

The criticality at the localization-delocalization transi- 
tion (LDT) is believed to be universal and entirely deter- 
mined by the generic symmetries of the systems (with 
certain special exceptions).^ One manifestation of this 
universality is a common scaling exponent v per generic 
symmetry class, describing the diverging behavior of the 
localization length^ at the mobility edge: <^ ~ {Ep-Ep)~^. 
This in turn leads to a temperature-scaling of the trans- 
port coefficients at LDT!^ 

, T) = F ({Ef - £^)(T/To)-^) , x = p/2v, (1) 

when combined with the concept of temperature- 
induced effective size introduced by Thouless,^ and with 
the single-parameter scaling hypothesis.^ ^ In Eq. ([l]), F 
is a system dependent function. To is a reference temper- 
ature, and p is the dynamical exponent for dissipation. 
Besides these universal scaling laws, there are other in- 
teresting universal aspects of the LDT, such as the exis- 
tence of a single critical point (as opposed to a line of 
critical points) for the unitary symmetry class, universal 
values of the critical transport coefficients and universal 
renormalization flow-diagrams of the transport coeffi- 
cients with the temperature or with the system-size. 

Such universal characteristics at the plateau-plateau 
(PPT) and plateau-insulator (PIT) transitions in the Inte- 
ger Quantum Hall Effect (IQFi E) ha ve preoccupied the 
experimental^ and theoreticaPl^^n^ condensed matter 
communities for decades, resulting in some of the best 
experimental data and computer simulations available 
for a quantum transition. With the discovery of Topo- 
logical Insulators (TI),^^^^ the principles of universality 
will receive renewed scrutiny. The TIs have bulk ex- 
tended states even in the presence of strong disorder,'^ 
hence they are expected to display sharp LDTs. As such, 
the transport measurements at the transitions could be 
as clean and revealing as the ones in IQHE. The TIs can 
fall in different symmetry classes^^^^ and, even within 



the same symmetry class, the topol ogica l materials can 
be very different from one anotherj^^Ml ^]^^g providing 
the perfect laboratory to test the principles of universal- 
ity. The computer simulations have already began this 
process (note that these are all zero-temperature finite- 
size scaling simulations) One question that received 
great attention from these works is if the scaling expo- 
nents of the symplectic models at the metal-to-normal 
insulator and at the metal-to-topological insulator are 
the same. So far, the answer seems to be affirmative. 

Although there is a substantial experimental transport 
data for TIs, only recently the focused was fully tuned on 
the topological-to-nontopological phase transition.EHS] 
And even these experiments need to be further refined 
for the quantum criticality to be revealed. Traditionally, 
the experiment was always ahead of the theoretical sim- 
ulations in this domain (see Ref . 35 for a discussion), but 
we strongly believe that this state of affairs will be soon 
reversed by the adoption of the non-commutative geom- 
etry program for aperiodic solids, initiated by Bellissard 
et al in the 90's.^^ For example, this natural formal- 
ism for treating disordered solids under magnetic fields 
enabled us to develop extremely accurate, stable and ef- 
ficient simulations of the zero- and finite-temperature 
transport coefficients,*^^*^^^ to a point where qualita- 
tive and quantitative agreement between experiment 
and simulation was possible for PIT in IQHE.^^ In this 
article, we announce several predictions based on the 
noncommutative Kubo formula on the quantum criti- 
cality at the Chern-to-normal insulator transition. On 
a broader note, we want to announce that this type of 
simulations reached a level where they provide qualita- 
tive and quantitative guidance for the experiments on 
quantum criticality at the topological-to-nontopological 
phase transitions. 

The search for possible Chern Insulator (CI) materials 
have intensified and several theoretical studies have al- 
ready singled out possible CI candidatesPE! This and 
the fact that the CIs and the IQHE liquids belong to the 
same unitary symmetry class motivated us to focus ex- 
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clusively on CIs in this study. For the spin-up sector of 
the Kane-Mele (KM) model^^ with strong disorder, we 
were able to converge {i.e. eliminate any finite-size ef- 
fects) the finite-temperature conductivity a and resistiv- 
ity p = tensors, at temperatures low enough to enter 
the quantum critical regime at the transition between a 
CI and a normal insulator. We compare these results 
with the known facts for the PIT in IQHE.^^ 1^ 22 24 82 85 

We find that, like at PIT in IQHE, the graphs of p^x as 
function of electron density, recorded at different tem- 
peratures, intersect at one single critical point, and they 
collapse into a single curve after a single-parameter 
rescaling. The scaling exponent k is in good agree- 
ment with what one would predict by using the univer- 
sally accepted value of the finite-size scaling exponent 
V = 2.58 + 0.03,^^^^ and with p fixed like in our sim- 
ulations (p = 1). We clearly see the expected renormal- 
ization flow of o with the temperature, but the separatix 
is not a semicircle like at PIT in IQHE, but rather an el- 
lipse. At the critical point, we find o^y ^ and, quite 
interestingly, a^x ~ j rather than | j (to be precise, the 
numerics place Oxy between 0.5 - 0.6 j at PIT^^^^^ ^ 
main surprise was, however, the absence of the Quan- 
tized Hall Insulator phase, characterized by cr = but 

Pxy = ^• 

To probe the broader applicability of our conclusions, 
we have repeated the computations for the spin-up sec- 
tor of the Berne vig-Hughes-Zhang (BHZ) model .^^ Al- 
though in this case the calculations are not converged 
well enough for accurate quantitative predictions, they 
already support the qualitative conclusions we reached 
for the KM model: existence of a single critical point 
and absence of the Quantized Hall Insulator pahse. The 
critical values of the transport coefficients remain the 
same. Unfortunately, the data for BHZ model display 
a poorer scaling, indicating that lower temperatures are 
needed. As such, the scaling exponent for this second 
set of calculations is inconclusive. We plan to re-examine 
the issue in future simulations with larger system-sizes, 
which will enable us to further reduce the temperature 
without introducing size effects. 



II. THE NONCOMMUTATIVE KUBO FORMULA 

The key for our simulations is the noncommutative 
Kubo formulaPS 

a,y(£f,T) = -r(feH](l/T + Xh)"'[X;,Ofd(H)]), (2) 

where T represents the trace over volume, H is the dis- 
ordered Hamiltonian, x is the position operator, t is the 
relaxation time for dissipation, Ofd is the Fermi-Dirac 
distribution for a given T and Ef , and Xh is the Liouvil- 
lian acting on the observables (i.e., operators): 

£hA = z[H,A]. (3) 



In real condensed matter systems: t - T ^ , where p is 
the dynamical exponent appearing in Eq. Q. 

As usual, the conductivity tensor cry^ of Eq. ^ pro- 
vides the link between the charge current density and 
the homogeneous static electric field E that drives the 
system: 

d 

}j = Y^aj,E,. (4) 

k=l 

Here, / is the time average of the expected value of the 
current density operator. The state of the system evolves 
according to the dynamics induced by H - eEx, and 
by random (in time), instantaneous electron scattering 
events. Realistic electron-electron and electron-phonon 
scattering matrices can be considered, however, Eq. ^ 
assumes the relaxation time approximation. 

Extended discussions of the physical assumptions, of 
the derivation and of the noncommutative formalism, 
together with convergence tests and applications to the 
disordered Hofstadter model can be found in Ref . [68l 
Another complete discussion of the Kubo formalism was 
given in this very recent journal (see Ref. |69|, together 
with applications to a disordered model of a quantum 
spin-Hall insulator. For these reasons, here we only dis- 
cuss the significance of Eq. (|2| and its relation with other 
works. Eq. (|2| represents the thermodynamic limit of 
the formal Kubo formulas written in terms of the eigen- 
functions and eigenvalues of the equilibrium Hamilto- 
nian, found in the classical solid state textbooks [see 
for example Eq. (3.385) and its finite-temperature ver- 
sion in Ref. ^97]. Finite-volume versions of the Kubo 
formula for aperiodic systems can be also derived us- 
ing the traditional quantum master equation.^^ In this 
type of analyses, the Kubo formula is formulated as the 
infinite-volume limit of the finite-volume expressions. 
The existence of this limit for aperiodic systems is vir- 
tually impossible to prove without the noncommutative 
formalism. What is special about the latter is that it en- 
ables one to work directly in the thermodynamic limit. 
For example, the trace per volume is rewritten using 
the Birkhoff ergodic theorem, and the Liouvillian can be 
defined as a normal operator on a well defined Hilbert 
space (hence a spectral decomposition exists for it). 

At the practical level, an explicit Kubo formula for 
infinite volume is important because it enables one to 
analyze how fast are the various finite-volume approx- 
imations converging in the thermodynamic limit. In 
particular, it enabled us to develop a canonical finite- 
volume approximation that converges exponentially fast 
in the thermodynamic limit, a fact that was established 
with mathematical rigor.^^ This canonical finite-volume 
approximation and its numerical implementation have 
been extensively discussed in Refs. 681 and 1691 Here we 
closely follow these two references. 
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FIG. 1. (Color online) (a) The spectrum of the disordered Chern 
insulator is localized ever3rwhere except for two thin energy re- 
gions, (b) The expected flow of the transport coefficients when 
the temperature is lowered to zero. Different flow lines corre- 
spond to different Fermi levels. The dashed line represents the 
separatrix, which is expected to satisfy a semicircle law. 



III. QUANTUM CRITICALITY IN THE SPIN-UP 
SECTOR OF THE KANE-MELE MODEL 

In the absence of Rashba interaction, the spin sectors 
decouples and the KM model reduces to a disconnected 
sum of two Chern insulators. The Hamiltonian in the 
spin-up sector (tuned in the middle of the CI topological 
phase) takes the form: 

Ho = Y, I^X^I + Yj ^n(\n){fn\ - (5) 

(nm) {{nm)) 

Here, t„ represents the isospin of the site n of the hon- 
eycomb lattice, and ()/(()) symbolize first/second near- 
est neighbors. We add the random potential V^j = 
W Yju to Ho, where the co's are independent ran- 

dom variables uniformly distributed in [-^, |]. We fix 
W = 4 (= 2 X the clean insulating gap), in order to achieve 
the strong disorder regime where the insulating (spec- 
tral) gap is closed and only a mobility gap remains. The 
simulations are performed on a finite-size lattice con- 
taining 80 X 80 unit cells. By repeating the simulations 
with different lattice sizes, we concluded that the effects 
due to the finite-size of the simulation box are practically 
negligible for the temperatures considered in the present 
study (the effective Thouless length^ is smaller than the 
simulation box). 

In the above conditions, all the quantum states of 
the model are localized, except for the states in two 
narrow energy regions separating the CI from the nor- 
mal insulator [see Fig. [T]ja)]. The existence of such de- 
localized quantum states can be demonstrated^^ with 
mathematical rigor using the theory of noncommutative 
Chern number,^^ while numerically, it has been demon- 
strated using rec ursive Gr een's function and transf er raa- 
trix calculations,SBl9M level statistics analysis,^^!^^ 
simulations of the edge currents and computations of 
the edge conductance. ^^^"^^^ Near the transitions, the 
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FIG. 2. (Color online) o as function of £f at /cT = 1/t = 0.01, 
0.02, 0.03, 0.04, 0.06, 0.08, simulated on an 80 x 80 lattice. An 
average over many disorder configurations was considered (as 
many as 67 for T = 0.01 and 24 for T = 0.08). The marks posted 
on the graph give a at the critical point. The inset shows the 
flow of our data in the (cr^y, Oxx) plane as T — > 0. The dotted 
line suggests the separatrix for the flow. 



field-theoretic arguments developed by Pruisken and 
collaborators^^^^^ for IQHE predict the T-driven flow- 
diagram shown in Fig. [ijb), which was observed and 
confirmed in IQHE by experiment.^^ 

The simulated cr is reported in Fig. |2] as function of 
Fermi level, for kT = 1/t = 0.01, 0.02, 0.03, 0.04, 0.06 
and 0.08 (hence we fix p = 1 in our simulations). In this 
figure we can see an energy region where, especially for 
the lower temperatures, Oyx takes the quantized value of 
-1, indicating that the system is in the CI phase. When 
moving away from this energy region, Oyx starts to con- 
verge towards 0, indicating that the system enters the 
normal insulator phase. The Oyx curves computed at 
different temperatures intersect each other at practically 
one point. Examining the data for the direct conduc- 
tivity, we see Oxx decreasing as T ^ for most part of 
the energy spectrum, a hallmark of the insulating phase, 
with the exception of two distinct energy regions where 
Oxx increases as temperature decreases. These energy 
regions appear exactly where Oyx switches between its 
quantized values and, as such, they must harbor ex- 
tended quantum states. The energy regions where Oxx 
increases as T ^ appear to become smaller and smaller 
as T is lowered, and that the maximum value of Oxx sta- 
bilizes at a finite value (as opposed to diverging to in- 
finity). An important question is if these regions reduce 
to a point as T ^ 0. A more refined analysis based on 
Fig. [S] shows that this is, indeed, the case, and gives the 
critical Fermi levels £^ ^ +1.02. The values of a at are 

marked in Fig. |2] and are crj^ ^ j and crj^ ^ ^j. These 
critical values for the conductance are reproduced when 
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FIG. 3. (Color online) (a) pxx (and p^y in the inset) as function of 
Ef at kT = 0.01, 0.02, 0.03, 0.04, 0.06 and 0.08. (b) p^^ as function 
of temperature for various £f values. The arrow indicates the 
transition from CI to the normal insulator. 



FIG. 4. (Color online) pxx at different temperatures collapses 
onto a single curve after the single-parameter rescaling: £f 
El + {Ep - £y(T/To)-^ with E^ = 1.017, Tq = 0.04 and k = 0.21. 



the simulations are repeated for the spin-up sector of the 
BHZ model. Furthermore, the PIT is known to exhibit 
a universal critical conductance,^^ though w^ith a differ- 
ent value, so we are led to conjecture that the Chern- 
to-normal insulator transition also exhibits a universal 
conductance, with the universal values stated above. 

When the data points from the main plot in Fig. |2] 
are placed in the (cr^^-y, cTjj) plane, we obtain the flow- 
diagram shown in the inset of Fig. 2] The separatrix for 
this flow, marked with the dottealine, has the shape 
of a semiellipse with the semiaxes and j. At PIT 
in the IQHE, the separatrix strictly obeys the semicircle 
law: al^ + (a^y - 03e^/h)^ = {O.be^lhf}^ Let us point out 
that the values Oxy > 0.5 occur when £f is located inside 
the old clean insulating gap, while the values Oxy < 0.5 
occur when £f is located in the high-density spectrum 
resulted from the localization of the clean energy bands. 
As such, the flow (as T ^ 0) starts from the inside (out- 
side) of the ellipse and moves towards the separatrix for 
Oxy > 0.5 {Oxy < 0.5), a markedly different behavior when 
compared with PIT. 

The existence of a unique critical point (as opposed to 
a line of critical points) at the Chern-to-normal transition 
can be established with great confidence by examining 
the resistivity tensor, plotted in Fig. [3] From the expres- 
sion Pxx = Oxxl{o^xx + ^jy)/ it follows that Pxx inside 
the CI phase, while pxx ^ inside the normal insu- 
lating phase, when T ^ 0. As a consequence, the pxx 
curves at different T must cross each other, very much 
like the Oyx curves do in Fig. |2] Since the plots are sym- 
metric relative to the zero energy, we can focus only on 
the positive energies. Fig. [3ja) strongly suggests that 
all Pxx curves cross each other at a single point, exactly 
how it was observed at PIT. To accurately pin-point this 
crossing point, we replotted pxx in Fig. |3lb), this time as 



function of temperature, for each positive Ef-value ap- 
pearing in Fig.|3ta). The flow of pxx with T curves down 
at lower E/s and it curves up at higher E/s. There 
is a clear separatrix between these two distinct tenden- 
cies, very much like in the experimental data on the PI 
transition in Ref . IllOi or on the metal-insulator transi- 
tion in Ref. Illll This enables us to accurately identify 
the critical point: E^ = +1.017, and then to determine 
the value of the conductivity tensor at the critical point: 
a'^y = -0.53 X f and aj^ = 1.09 x f . 

The inset in Fig.|3];a) reports pyx = Oxyl{ol^ + oly) as 
function of £f , which decreases from 1 to almost with 
the same rate as Oyx. As one can see, there is absolutely 
no tendency for pxy to stay quantized at ^ through the 
transition or further into the normal insulating phase. 
Such quantization of pxy would have been incompatible 
with the critical values of o determined above. It is 
also known that a quantized pxy is equivalent with the 
semicircle law,!^ but the separatrix shown in the inset 
of Fig. [2] has an elliptical shape quite different from a 
semicircle. We want to point out that these facts are 
also true for the simulations with the BHZ model. As 
such, we can conclude with great confidence that the 
Quantized Hall Insulator phase is absent for this system. 
This is in striking contrast with the PIT in IQHE, for 
which we did observed the Quantized Hall Insulating 
phase,^^ using same type of calculations. 

We now zoom into the region around E^ and start the 
scaling analysis. Since the scaling occur in the asymp- 
totic limit T ^ 0, we keep for this analysis only the low- 
est four temperatures. As shown in Fig.lij the pj^^-curves 
collapse almost perfectly on top of eacn other after the 
energy axis is rescaled as: Ef ^ £^ + (Ef - £^)(fcT/fcTo)~^ 
(fcTo = 0.04). The best overlap of the rescaled curves 
is obtained for k = 0.21 + 0.01, a value that is in good 
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FIG. 5. (Color online) The simulated o for BHZ model as func- 
tion of Ef at /cT = 1/t = 0.01, 0.02, 0.04, 0.06, 0.08, simulated on 
an 80 X 80 lattice. The inset shows the T-driven renormalization 
flow of o in the (cJxy, Oxx) plane. 



agreement with k = 0.194 + 0.002 obtained from the ex- 
pression K = p/(2v) with the universally accepted value 
V = 2.58 ± 0.03, and p = 1 like in our simulations. 



IV. CRITICAL REGIME IN THE SPIN-UP SECTOR OF 
THE BERNEVIG-HUGHES-ZHANG MODEL 

This section reports the simulations based on the non- 
commutative Kubo formula, for the Chern Insulator (CI) 
corresponding to the spin-up sector of the Bernevig- 
Hughes-Zhang (BHZ) model (tuned in the middle of 
the topological phase) 

h{k) = Ox sin kx + Oy sin ky + 2(1 + cos kx + cos fcy)^^, (6) 

where the cr's represent the Pauli's matrices. This model 
can be represented on a 2-dimensional square-lattice 
with two orbitals per site. The lattice sites are indexed by 
n and the orbitals by a. In this real-space representation, 
we add the random potential Vo) = W X,^,^ 
where c^^^ creates an electron in state a at site n, and 
co's are independent random variables uniformly dis- 
tributed in [-1/2, 1/2]. We fixed W = 5{=2.5x the clean 
insulating gap), to achieve the strong disorder regime 
where the insulating gap is closed and only a mobility 
gap remains. As in the previous section, the lattice-size 
was taken to be 80 x 80 unit cells. An average over 
many disorder configurations was considered for dif- 
ferent temperatures, specifically: 68 configurations for 
kT = 0.01, 67 for kT = 0.02 and 23 for kT = 0.04, 0.06 and 
0.08. 

Fig. |5] reports the simulated conductivity tensor a as 
function of Fermi level Ep, at various temperatures T. 
The data display the same general features as seen in 
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FIG. 6. (Color online) (a) pxx for BHZ model as function of Ep 
at kT = 0.01, 0.02, 0.04, 0.06 and 0.08. The inset reports p^y. 
(b) pxx as function of temperature for various Ep values. The 
arrow indicates the transition from CI to the normal insulator. 



Fig. |2] In particular, Oxy takes quantized values in cer- 
tain energy regions and, in the energy region where Oxy 
shifts between the quantized values, Oxx becomes inde- 
pendent of temperature, signaling the presence of ex- 
tended states. A more refined analysis based on Fig. [6] 
reveals again a single critical point, located at the critical 
Fermi level ^ 0.94. The critical cr values at £^ were 

found again to be crj^ ^ j and cTjy ~ | f (the exact values 
are marked on the graph). The inset of Fig. |5] shows the 
T-driven flow of a in the {Oxy^ Oxx) plane. While it gener- 
ally resembles the flow in Fig.|2j a closer analysis reveals 
that the calculations are not well converged on the nor- 
mal insulator side (near Oxy = 0). Lower temperatures 
will be definitely needed to obtain the correct separatrix 
(which we believe it will look very similar to that for KM 
model). 

Fig. [6] reports the simulated resistivity tensor p (a) as 
function of £f , at various temperatures, and (b) as func- 
tion of fcT, at various Fermi levels. The flow with T 
in Fig. [6tb) is used to accurately determined the critical 
point, whose coordinates were given above. Apart from 
some structure in the curves, there are no major differ- 
ences when Fig. |6] is compared with Fig. |3] The inset 
reports pxy as function of Ep and, here too, there is no 
trace of the Quantized Hall Insulator phase. 

Fig. [tI reports the scaling analysis for the BHZ model. 
Only trie lowest three temperatures have been consid- 
ered. The best overlap of the rescaled curves is obtained 
for K = 0.14+0.01, but as already anticipated, the overlap 
is quite poor. As such, we conclude that this numerical 
value is inconclusive. 

When searching for a reason for the poorer conver- 
gence of the results in the BHZmodeL we found that, 
in all our previous 

simulationsEEzEIm 

Uheband energy 

states in the BHZ model localize much slower than in 
the KM model. The only major qualitative difference 
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FIG. 7. (Color online) Single-parameter rescaling for BHZ 
model: £f ^ £^ + (£f - El){T/To)-'', with £^ = 0.94, To = 0.01 
and K = 0.14. 



between KM and BHZ models is that the latter has one 
split-Dirac point, located at A: = 0, while the former has 
two split-Dirac points located at fc 9^ 0. This may indeed 
alter the physics of impurity scattering processes, re- 
sulting in the distinct behaviors that we observed in our 
present study What is certain is that lower temperatures 
are needed to fully enter the critical regime in the BHZ 
mode, something that we defer to future investigations. 



V. CONCLUSIONS 

In conclusion, the simulations based on the non- 
commutative Kubo formula and a recently developed 
numerical implementation enabled us to converge the 
transport coefficients at temperatures low enough to en- 
ter the quantum critical regime at the Chern-to-normal 
insulator transition, at least for the Kane-Mele model. 
When compared with the available experimental facts 
and our previous simulations for PIT in IQHE, the re- 
sults on the two strongly disordered Chern insulator 
models show similarities but also important differences. 
The similarities include: the existence of a single criti- 
cal point and the single-parameter scaling behavior; The 
KM model, for which the full critical regime seemed to 
be reached by our calculations, displays a scaling expo- 
nent consistent with the universally accepted value for 
the unitary class. Among the dissimilarities were the ab- 
sence of the Quantized Hall Insulator phase, a universal 

critical value of crj^ ^ j instead oi a^^^ ^ and the 
violation of the semi-circle law. 
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